Dynamically dominant excitations of string solutions in the spin-1/2 
antiferromagnetic Heisenberg chain in magnetic fields 
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Using Bethe-ansatz solutions, we uncover a well-defined continuum in dynamical structure factor 
5*^ (fc, Lo) of the spin-1/2 antiferromagnetic Heisenberg chain in magnetic fields. It comes from string 
solutions which continuously connect the mode of the lowest-energy excitations in the zero-field limit 
and that of bound states of overturned spins from the ferromagnetic state near the saturation field. 
We confirm the relevance to real materials through comparisons with experimental results. 
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The spin-1/2 antiferromagnetic Heisenberg chain ex- 
hibits intriguing quantum many-body effects, associated 
with modern concepts in condensed-matter physics, such 
as spin liquids, quantum criticality and fractionalization. 
Also, this system is an excellent platform to make pre- 
cise comparisons between experiments and theories: Var- 
ious intriguing features predicted by exact solutions [ij 
can be confirmed by accurate experiments on quasi-one- 
dimensional materials. Actually, it is established that 
dynamical properties in the absence of magnetic field 
are characterized by quasiparticles called spinous [1, 0| 
through precise comparisons between theoretical predic- 
tions and experimental results 0, H, 11] ■ 

In magnetic fields, dynamical properties are more com- 
plicated. Some basic features of dominant excitation 
spectra can be understood by modifying the 2-spinon 
continuum in zero field accordingto the strength of mag- 
netic fields in the k-uj plane 0) fi @| • The distributions 
of spectral weights in S ^{k,uj) and S^^{k,uj) are ef- 
fectively expressed by low-order par ticle-hole excitations 
in Bethe-ansatz solutions [3, [lol. Illl| : The dominant ex- 
citations in S ^{k,Lu) and S^^{k,uj) are known as 2- 
psinon {2ip) excitations and psinon-antipsinon (ipip*) ex- 
citations, respectively [lol 11 1. Their properties have 
been inve stig ated in detail by using Bethe-ansatz solu- 



tions 



3, ETIil, H Q E S [13, [3 • 

{k,uj) 



As for {kjUj), the situation is much more compli- 
cated. Except the low-energy modes near momentum 
fc=0 [i, H, [S IB, [2^ and k=Tr d 0, E , behav- 
iors of have not been clarified: In Ref. 0, the 
continuum of 'llJ^p* excitations was predicted to persist 
in the thermodynamic limit based on the classification 
by the Bethe formalism. However, numerical results in- 
dicated that the spectral weight in the tpip* continuum 
is rather small except near the edges, and there exists 
a large fraction of spectral weights above the continuum 

In this Letter, mainly focusing attention on behaviors 
of in magnetic fields, we identify excitations 

having large spectral weights to clarify overall dynamical 
features of the Heisenberg chain in magnetic fields. 

We consider the spin-1/2 antiferromagnetic Heisenberg 



chain with L sites, M down spins {M<L/2) and periodic 
boundary conditions. The Haniiltonian is defined as 



H — J ^ Sx ■ Sx+i — HS^ 



(1) 



where Sx is the spin-1/2 operator at site x, and J>0. 
Magnetic field H at S"" / L=l/2-MjL in the thermo- 



dynamic limit is obtained in Ref. [21|. In the Bethe 
ansatz [ij, the wave function, energy and momen- 
tum of an eigenstate are expressed by a set of rapidi- 
ties {Aj} which is obtained from the Bethe equation: 
L arctan( Aj ) =7r/j - 
of {Ij} is given. 



-X)fiia'rctan[(Aj— A;)/2], once a set 
Here, Ij {j=l^M) are called Bethe 
quantum numbers, which are integers (half-odd integers) 
within |/j |<(L-M-l)/2for odd (even) L-M as in Fig.[T] 
(a-e) for solutions of real {Aj}. Distributions of {Ij} are 
somewhat analogous to momentum distributions of spin- 
less fermions. As in Fig. [T](b), a hole (particle) created 
inside (outside) the {Ij} of the ground state is called psi- 
non (antipsinon) and denoted by -tp (ip*) [13, [Hi- 
lt is also known that there are solutions with complex 
rapidities [l|. Later, we will consider solutions with an 
n-string (n>2), i.e. a set of n rapidities expressed as 
Aj=A+i{n+l—2j)+irij for j=l^n, where A is real, and 



rij^Oie-"^) with oO P,[3,|2i|23|. We take real A^ for 
i>n+l, and denote the n-string solutions by Sn. These 
solutions are specified by two sets of {Ij} [22]: One is for 
the real rapidities, and the other is for the n-string, which 
we denote by {/j} and {I]}, respectively. For Sn (n>2), 
{I')}^I\, \Ii\<L/2-M as in Fig. [1] (f-2). We denote tp 
and V'* of {Ij} by ips and ip*, respectively, and regard 
{Ij} without ips or ip* also as a part of those with ipsi'l- 
We investigate behaviors of dynamical structure fac- 
tors defined as 5""(fc, w)=X;i Af""(fc, e^)(5(w - e^) for 
a=— , -I- and z. Here, M°^°'(k,e\) is the transition rate 
between the ground state |G.S.) in a magnetic field and 
an excited state |fc, e\) with excitation energy e\. and mo- 
mentum A:, defined as M""(fc,e*fe)=|(fc,eyS'^|G.S.)p. We 
calculated M°"^{k, ej) of the Heisenberg chain, following 



Refs. [il [il [13, [3, [53, [13 where M""(fc,e^) is ex- 
pressed in a determinant form whose matrix elements are 



/, -8-7 -6-5 -4-3 -2-1 012345678 
,/ I I I I I I I I I I I I I I I I I 
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FIG. 1: Distributions of Bethe quantum numbers {7j }. Filled 
symbols denote {Ij} of the ground state of Af=8 and typical 
excited states in dynamical structure factors for 0<fc<7r in 
I/=24. Blue open square and red solid diamond denote psinon 
(■0) and antipsinon JUJ], which are located at points on 

thick and thin lines, respectively, (a) Ground state, (b) tpip* 
for S'^^(fc,tj). (c) 2V'for5'"+(fe,tj). (d) 2?/)* for the continuum 
near k—TT in S^^ {k, cj). (e) 1?/)* for the low-energy mode near 
k=0 in S'^^ {k,uj). (f-1) lip of {7j} for 2-string solutions in 
S+"(fc,w). (f-2) Ji for the n-string of Sn (n>2) in S+-{k,u}). 



expressed in terms of rapidities [24]. As for string solu- 
tions with an n-string (n>2), we calculated M""(fc,e5.) 
after transforming the matrices so that sin gula rities in 
determinants can be cancelled as in Refs. jlil . [23j . In 
(fc=0, oj), we took into account the resonant mode of 
the field-induced magnetization [^[13, El, 3 2^- In this 
Letter, we show results on dynamical structure factors in 
units of 1/J for 0<fc<7r, noting 5""(fc, t^)=S'""(-fc, w), 
a=— , + and z. 

We calculated the spectral weight in L=320, using up 
to 2'ip2ip* excitations, and evaluated the contributions 
to sum rules '25'|. The results are shown by open sym- 
bols in Fig. [2] (a) and (b). The spectral weight of up to 
2'02V'* excitations in S ^{k,uj) satisfies more than 90% 
of the sum rule in L=320 as shown by open diamonds in 
Fig. [11(b). The spectral weight in 5*^ {k,uj) severely de- 
creases in the small regime as shown by open circles in 
Fig. [1(a). Figure [2] (e) shows S+-{k,uj) in the smah S"" 
regime. Although (fc, to) should be continuously con- 
nected to that in zero field (Fig. [2| (d)) as the magnetic 
field decreases, most of the spectral weight in {k,uj) 
is missing. In particular, the large intensity near the des 
Cloizeaux-Pearson mode (dCP) [2^ (the lower edge of 
the continuum in zero field) is almost lost. This implies 
that there are other important excitations for (fc, uj). 

To explain the origin of the missing spectral weight, we 
consider S2 and S3 defined above. Among many states 
in S2 and S3, we found that S2 with V'sV's ^.nd S3 with 
2'0s have large weights in {k,uj). By taking into ac- 
count these string solutions, ik,uj) in the small 
regime almost recovers the missing weight near the dCP 
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FIG. 2: (a,b) Contributions to sum rules from solutions of up 
to 2ip2il)* , string solutions with a 2-string and those with a 3- 
string, denoted by R, S2 and S3, respectively, in L=32Q for (a) 
5'+-(fc,w), (b) S- + ik,u}) and S^^(it,tj). (c) Contributions of 
dynamically dominant excitations in L=2240. (d) S~'^{k,uj) 
in zero field. (e,f) S'*'~{k,iu) at S^/L=0.025 without and with 
the contributions of S2 and S3. In (d-f), the data obtained in 
I/=320 are broadened in a Lorentzian form with full-width at 
half maximum (FWHM) 0.08 J. For S2 and S3 in (a,b,f), we 
took the O(L^) states described in the text. 



mode and that of the continuum near k—ir as shown in 
Fig.[2|(f). The intensity near the dCP mode is mainly due 
to the S2, and that of the continuum near k=TT is mainly 
due to the S3. The contributions from these string solu- 
tions increase as the magnetic field decreases as shown by 
solid symbols in Fig.[H(a). By taking them into account, 
more than 80% of the total spectral weight in 5"' {k,uj) 
is explained in a wide range of in L=320 (27j . 

As for n-string solutions (n>2) in S^^(k,uj), we found 
that a large spectral weight is carried by S2 with 2-08 in 
low fields. As shown by solid squares in Fig. [D (b), the 
decrease of S^^{k,uj) of up to 2ip2ip* excitations in the 
small regime (open squares in Fig. [D (b)) is almost 
compensated for by the S2, and more than 90% of the 
total spectral weight in 5^^(fc,a;) is explained in a wide 
range of in i=320 [27]. These string solutions are 
relevant to the continuum near k=TT in low fields. 

Figure [3| (a) shows behaviors of dynamical structure 
factors in magnetic fields in L=320, where we took 
into account excitations of up to 2ip2ip* and the above- 
described string solutions. Hereafter, we use these ex- 
citations to calculate dynamical structure factors. In 
(kjOj), there are three sets of dominant continua as 
shown in the top row of Fig. [3] (a). The high-energy con- 
tinuum (a;/J>2), which is mainly due to S2, goes up to 
higher energies as the magnetic field increases, separated 
from the low-energy continua. Although this continuum 
has been observed in the high-field regime by numeri- 
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FIG. 3: (a)S'+'(fc,cj)/2, S" + (fc,a;)/2, S^^{k,iu) and ^""(fc, a;) 
(from above) at S'=/L=l/8, 1/4 and 3/8 (from the left) in 
L—320. The data are broadened in a Lorentzian form with 
FWHM=0.08J. (b) Same as (a) but those of dynamically 
dominant excitations in L—2240 with the product ansatz. 



cal calculations [3, Hi > the relation to the S2 has not 
been clarified. Near the saturation field, these string solu- 
tions reduce to bound states of two overturned spins from 
the ferromagnetic state [H with the excitation energy 
given by a;(fc)=-^(3— cos A;)-|-iJ. As for S ^{k,uj) and 
S^^{k, Lo), there is basically one continuum for each corre- 
lation as shown in the second and third rows of Fig. [3] (a), 
respectively. The additional high-energy continuum at 
S^/L— 1/8 in S''^{k,uj) near k=TT mainly comes from S2. 

From many excitations considered above, we extract 
dynamically dominant excitations of 0{L^) states which 
characterize the behaviors of Fig. [3] (a). Then, the 
product ansatz [111: S""{k,Lo)=M""{k,Lu)D{k,uj) with 



Dik,el)=2/{e 



e]r^) can be applied. Here, in energy 
regions where there is more than one sequence of states, 
we took the one with the largest spectral weight. Here- 
after, we denote S'""(A:,a;) of the dynamically dominant 
excitations as 'S'j^(fc,cj) for a=— , + and z. Their con- 
tributions to sum rules in L=2240 are shown in Fig.[2](c). 

The top row of Fig.E](b) shows Sj~^{k,u}). The high- 
energy continuum near w/ J>2 is due to S2 with IV's as in 
Fig.[l](f-1), which can be regarded as a part of the above- 
mentioned S2 with ipsijjs- The low-energy continuum 
near fc=7r comes from 2-0* excitations where two ip*s are 
located on opposite sides of the filled region of {Ij} from 
each other as in Fig. [1] (d). The low-energy mode near 
k=0 is due to excitations with one ^* in the right empty 

■ [^e). This mode has 
1( 



region of {Ij} for 0<fc<7r as in Fij 
been mentioned in the literature ^ 

The dynamically dominant excitations for S ^{k,uj) 
and S^^{k,U!) are known as 2ip and ^|njJ* excitations as in 
Fig. [11(c) and (b), respectively 13, The results are 
shown in the second and third rows of Fig. [3] (b) . 

Noting that 2ip, ipip* and 2ip* excitations are dynam- 
ically dominant in S ^{k,uj), S^^{k,uj) and {k,uj) 
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FIG. 4: Comparison with experimental results on 
CuCl2-2N(C5D5). (a) Solid lines in upper two panels are the 
present results of S^^{k,ijj) at fc=7r/2 in L=320. The data 
are broadened in a gaussian form with FWHM~0.17J. Open 
symbols are experimental results in Ref. 28]. The height is 
rescaled after subtracting the background (=180 counts). The 
lowest panel shows lineshapes of S^^^{k,ij) (a=— , + and z) 
at fc=7r/2 in L=2240. (b) SZn-,{k,i^) in L=2240. (c) Same as 
the middle panel of (a) except that the solid line is obtained 
from solutions with real rapidities only. The arrow indicates 
the high-energy peak observed in the experiment [28| . 



and that these excitations have 5*^=+!, and —1, re- 
spectively, we can naturally assign S'^=+l/2 and —1/2 
to Tp and ip* , respectively. Also, noting that excitations 
of S2 with V'sV's and those with 2tps have large weights 
in 5*^ (kjOj) and S^^{k,uj) and that these excitations 
have 5'^=— 1 and 0, respectively, we can naturally inter- 
pret the quasiparticle representing the 2-string (Fig.[l](f- 
2)) as a bound state of two tp*s which carries 5"^=— 1. 
This assignment is also applicable to 4-spinon states of 
S ^{k,u!) in zero field which have four spinous and a 2- 
string. The If/;* mode near k—0 in {k, lo) can also be 
regarded as the most dominant part of a 27/;* continuum 
where both are located in the right empty region of 
{Ij} (Fig. [T](e) as compared with (d)) for 0<fc<7r. But, 
since the spectral weight in this continuum is very small 
except the lip* mode, the tp* in this mode practically be- 
haves as a quasiparticle carrying 5*^=— 1, which reduces 
to a magnon carrying S^=—l above the saturation field. 

To confirm the relevance to real materials, we com- 
pare the present results with available experimental data 
on quasi-one-dimensional materials which can be effec- 
tively regarded as spin-1/2 antiferromagnetic Heisenberg 
chains. In inelastic neutron scattering experiments, a 
quantity proportional to S^^{k,uj) is observed. Here, we 
define S''''{k,uj) = [S-+{k,u)+S+-{k,uj)+4S''{k,uj)]/6. 
The fourth row of Fig. [3] shows the results of S'^'^{k,uj) 
and S'jQ„j(fc,w), where S^l^{k,uj) denotes S'^"'{k,uj) of 
the dynamically dominant excitations. 

In Fig. m (a), we compare the present results with ex- 
perimental results on CuCl2-2N(C5D5) (CPC) We 
broadened the numerical results of S^^^kjUj) in a gaus- 
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FIG. 5: Comparison with experimental results on 
Cu(C4H4N2)(N03)2. (a) SZra{k,Lj) at 5'7L=0.15 in 
L=2240. (b) Solid blue lines denote S'"{k,uj) in L=320 along 
green dotted lines in (a). The data are broadened in a gaus- 
sian form with standard deviation a"=0.1 meV. Open symbols 
are experimental results of intensity at H=8.7 T in Ref. ^. 
Here, we subtracted the background equal to that of H=0 H . 

sian form in zero field, and rescaled the experimental data 
in Ref. [28*1 after subtracting the background. Using the 
same broadening and rescaling parameters, we calculated 
S'^^{k, uj) at iJ«70 kOe and compared it with experimen- 
tal results, where we used g-factor g=2.Q8 and J=27.32 
K The positions and intensities of the two peaks 

at H=70 kOe observed in the experiment are reasonably 
reproduced by the present results as shown in the mid- 
dle panel of Fig. 2] (a). The high- and low-energy peaks 
are mainly due to S2 in ik,Lu) and 2'ip excitations in 
S ^(A;,a;), respectively, as shown in the lowest panels of 
Fig. m (a) and (b). For comparison, the results of solu- 
tions with real rapidities only are shown in Fig. |4] (c). 
Obviously, the spectral weight at the high-energy peak is 
missing, if n-string solutions (n>2) are neglected. 

Figure (b) shows comparisons with experimental re- 
sults on Cu(C4H4N2)(N03)2 (CuPzN) ^ along scans de- 
noted by green dotted lines in Fig. [5] (a). The present 
results reasonably agree with the experimental results. 
Figure [5] (a) is plotted so that it can be easily compared 
with Fig. 4 (a) in Ref. @ . A signature of the high-energy 
continuum originating from the 2-string solutions was ob- 
served near fc~0.37r and ti;~2 meV in the experiment 

Signatures of the 2-string solutions in an anisotropic 
two-dimensional system will be shown elsewhere (30| . 

In summary, we have investigated dynamical proper- 
ties of the spin-1/2 antiferromagnetic Heisenberg chain in 
magnetic fields, using Bethe-ansatz solutions. We found 
that string solutions with a 2-string have considerable 
spectral weight in (k,uj), and can be regarded as 
dynamically dominant excitations of (fc, lu). This in- 
dicates that not only and "0* but also the quasiparticle 
representing the 2-string plays an important role for the 
dynamical properties in magnetic fields. The continuum 
of these solutions continuously connects the dCP mode in 



the zero- field limit and that of bound states of overturned 
spins from the ferromagnetic state [l[ near the saturation 
field. Another finding is that 2ip* excitations are identi- 
fied as the dominant excitations of near fc=7r 
in the low-energy regime, which leads to a natural inter- 
pretation of tjj and ip* as quasiparticles in magnetic fields 
carrying S'^=-|-l/2 and —1/2, respectively. Comparisons 
with available experimental results reasonably support 
the relevance to real materials. For clear identification of 
dynamically dominant excitations, further experiments 
in higher fields at high energies are desired. Behaviors 
shown in this Letter are expected to be more or less true 
for general spin-1/2 XXZ chains in magnetic fields. 
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